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Abstract 



We consider a variant of the Georgi Glashow model in the BPS limit, augmented by a higher derivative 
Skyrme-like term, which is the simplest YMH model that can support monopole bound states. The 
spherically symetric solutions are studied with a combination of analytic and numerical techniques, 
^ ^1 which strongly suggest that the solutions converge to a finite energy configuration in the limit of infinite 

^ ' coupling of the Skyrme-like term. 

^ ■ 1 Introduction 
> 

Q>^ ■ The asymptotic analysis for the unit charge 't Hooft-Polyakov monopole [H [2] was carried out long ago 

I by Kirkman and Zachos [3] and by Gardner [3], and, has recently been elaborated by Forgacs et. al. [S] 

lO ' by providing a high precision numerical analysis of the problem. In contrast to the Prasad-Sommerfield 

I monopole |B] which is evaluated in closed form, the 't Hooft-Polyakov monopole can be evaluated only 

CSl ' numerically since it is a solution to the Georgi- Glashow model that exhibits a symmetry breaking Higgs 

self interaction potential. The existence of this numerically evaluated solution is underpinned by the purely 
analytic proof of existence given by Tyupkin, Fate'ev and Schwarz [7], but this does not shed any light on the 
behaviour of the solution as a function of the strength A of the Higgs potential term. Numerical studies of 
this A dependence, in [5] for the spherically symmetric case and in [3] for the axially symmetric case, reveal 
in particular that as A ^ oo the energy of the monopole asymptotes to a finite value. It is this behaviour 
■ on A found numerically that is underpinned by the analytic analysis of [3] . 

I It is our intention in this short note to supply an asymptotic analysis similar to that of [3], for the unit 

charge monopole of a variant of the Georgi-Glashow model characterised by the addition of a Skyrme like 
term in terms of the covariant detivatives of the Higgs field. The model is the SO (3) Higgs model 

n = l QlJ^f P + \D^r\' + \n\Di,rD,]cl^'\'^ > 4ng^ , (1) 

K giving the strength of the coupling of the Skyrme like term, and k having the dimension of L*. The lower 
bound on the right hand side is the usual monopole charge density 

Qi = ^eaa'be'''F,fDkC^'' (2) 

featuring a Skyrme like term in lieu of the Higgs symmetry breaking potential in the Georgi-Glashow model. 
This monopole, which was described in [lOj as a Skyrmed monopole, is a solution to a model that was 
distilled from a rather more involved set of models [HI [12] , the latter being designed to support monolpoles 
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with both mutually attracting and repelling phases. (The Skyrnied monopole of [10^ supports bound states 
which are axially symmetric for charges higher than 2, in contrast to Skyrmions of charges 3 and higher 
which instead exhibit Platonic symmetries |13|.) 

In contrast to the 't Hooft-Polyakov monopole for which there is an analytic proof of existence [7j, 
there are no such proofs for the monopoles of the various generalised models [T^ [TP] . This is because 
of the presence of higher order Skyrmc like terms, and the only known solutions are those constructed 
numerically. Here we will supply an asymptotic analysis analogous to that of [5], for the unit charged 
Skyrmcd monopole which is the simplest such example available. Due to the considerably more complex 
structure of the equations here, it is very difhcult to perform the promised asymptotic analysis using purely 
analytic method, and intcad we present a combination of both analytic and numerical analysis. (Similar 
techniques were used in [14J, in the context of the 5*0(3) gauged Skyrmion [TB] .) 

In the next two sections we present the asymptotic and the numerical analyses, respectively, followed by 
a brief summary of our result. 



2 Asymptotic analysis 

Since we are restricting to the charge- 1 monopole of ([T]), we impose the usual spherically symmetric Ansatz 

^aa'] ^ 1-W{r) , ^,,] ^ ^ ^ ^^^^ 

where the brackets [ab] imply antisymmetrisation, x'^ is the unit position vector and rj is the dimensionful 
VEV of the Higgs field. Imposing ([3]), the residual one dimensional static Hamiltonian is 

H = w''+ ^""^ + ^ h'' +w'h' + Kw'h'( 2h'' + ^] , (4) 

Zr^ 2 \ ^ J 

having rescaled r ^ rjr and k — » ryK so that both the rescaled radial variable r and the rescaled Skyrme 
coupling K in (|4]) are dimensionless. The corresponding equations for w and h are 

w" + -^—^ - wh^ = 2Kwh^ ( [h'f + — - 



r'^h" + 2rh' - 2w'^h 



Anw^h (^{h')^ + - 4k {w^h^h" + 2w^h{h'f + 2ww'h^h') 



(5) 



Due to their complex structure we are only able to extract from the equations some asymptotic informa- 
tion near the origin, at infinity, and for k oo. The asymptotic analysis will then be complemented by the 
numerical results. For r — > 0, we have w 1 and — > 0. A dominant balance analysis gives 

w = 1 + W2r^ + W4»'^ + 0{r^) , h = hir + h^r^ + O(r^) 

By induction we see that the asymptotic expansions contain only even or odd powers of r for w and h, 
respectively. For k = 0, the Prasad-Sommerfield solution yields W2 — — ^ and hi = i. For nonzero k, W2 
and hi have to be determined numerically. The coefficients for the next highest order are 



3 2 1 /-, . ,2n,4 , 2-AKhl 



and all other coefficients can be calculated recursively. 

For r ^ oo, we have w — > and h ^ \. Here the dominant balance analysis leads to an exponential 
fall-off for w and to 
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^h" + 2rh' = 
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for the leading term in h. We therefore have 

h = l-'^ + (a(r)e-^'') , w = P{r)e''' + O (tWc^^'') (6) 

Using again induction, we see that the asymptotic expansion at infinity for h contains only even powers of 
, whereas the asymptotic expansion for w contains the odd powers of . The coefficient functions in 
front of the exponential functions, starting with a(r) and /3(r), are polynomially bounded. For k — 0, we 
have q = I and /? = 2r. For nonzero k, P(r) satisfies the equation 

/3"(r) - 2/3' (r) + ^Pir) + - ^) P{r) = 2« f ^ " ^ + 4) /^M 

and therefore 



3 Dependance on 

To study the dependence of the energy 



E ^ H dr 







on K, we calculate 

dE fdH dw dH dw' dH dhdH dh' dH 



dn ./n \ dK Ok? dw dn dw' dn dh dn^ dh 

w^h^ 



dr 







^^h'l2{hr + ^]dr''^Es,>0 (7) 



Here we have used integration by parts, the equations for w and h, and the boundary conditions. We see 
that the energy increases with k. We also see that, if the energy is bounded as k — > oo (as suggested by the 
numerical results reported in the next section), Esk and therefore the product wh must vanish in this limit. 
Because of the boundary conditions, w cannot be zero for small r, whereas h cannot vanish for large r. This 
strongly suggests that h is zero in some interval (0,rm) and that w is zero in the interval (r^, oo). 
In the limit k oo we therefore expect the following equations to hold, 

w"ir) + &^^^0, %)=0 (0 (8) 
w{r) = , r'^h"{r) + 2rh' [r) = (r,„ < r < oo) (9) 

with boundary conditions 

w{Q) — 1 , w{rm) — , h(rm) — and h ^ 1 as r — > oo 
The solution of equation ^ is 

h{r) = (0 < r < r™) , h{r) = 1 - ^ (r„ < r < oo) 

r 

The solutions to equation ([5]) have been studied a long time ago [T7], but not to the same extent as other 
special functions. One property we can deduce immediately from equation ([8|) is that, because w" is negative 
for r KTm, w'irm) cannot be zero. 

Denoting the solution of Eqs. ([H]),® with the appropriate boundary condition at r — Vm by Woc^oo, 
the corresponding energy can be obtained easily : 

E^ir^) = r f-S + dr + T + h'A dr = E^ + ^ ■ (10) 
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Figure 1: The profile for the solution for several values of k. 



From the second term we see that Eao ^ oo for r„i — + and for oo; so Eao{rm) has a minimum. 

Solving Eq. ([5]) numerically for several values of r„j and computing the value -Eoo(^m)j we have determined 
the local minimum of Eao which we find to occur for r^.c ~ 2.0623. Furthermore, our result strongly suggests 
the relation -Eoo (''m,c) — Tm,c] we have no analytic proof for this to be an identity but it holds whithin our 
numerical accuracy i.e. 10^''. The configuration minimizing Eoo has u'^(O) « —0.656. 

The results of the numerical analysis reported in the next section will strongly confirm that 

lim (w(r), /i(r)) = (woo, ^oo) with r„i = r^.c- (11) 

K — >oo 

for the solution w, h of Eqs. ([5]). 



4 Numerical analysis 

We now discuss the numerical solutions of Eq. ([5]) for finite n. In the limit k = 0, the classical equations 
coincide with the equations of the BPS monopole which is known explicitly. To the best of our knowledge, 
no explicit solution exist for k > 0. We have solved equations (O completed with the boundary conditions, 

w{Q) = 1 , h{Q) = , w{oo) = , h{oo) = I (12) 

by using a numerical solver [15]. The PBS monopole gets smoothly deformed for k > 0. This is illustrated 
in Fig. [U where the profiles w, h of the BPS monopole solution are superposed with solutions corresponding 
to several positive values of k. For reasons explained in the previous section, we supplemented this figure 
with the profiles of the product wh. 

Several parameters characterizing the solutions, namely the classical energy E^ the energy of the Skyrme 
term Eg^ (see (O), the value of q (see (O), hi = h'{0) and W2 = w"{0)/2 are plotted as functions of k 
in Fig. [21 The characteristics of the PBS monopole are recovered in the limit k = 0, for instance E = 1, 
q = 1, h'{0) = 1/3, w"{0) = —1/6. The natural challenge is to construct numerically the solutions for large 
values of k and to confirm that they evolve according to the pattern discussed in the previous section. The 
quantities W2, hi, q, E extracted from our numerical solutions are reported on Fig. [21 The figure shows that 
they stay finite for k 1. In particular for k > 10^ we find q « 2.0623, w"{0) « —0.656, h{0) « suggesting 
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that, in the hmit k ^ oo, the solution approaches the configuration ([S]),® which minimizes the energy E^x^ . 
Our numerical results further indicate that the the product wh tends uniformly to the null function in the 
limit K — > oo. As anticipated in the previous section, we observe from our numerical solutions that, in the 
interior region, i.e. for r G [0,rm] we have h{r) ~ 0, while in the exterior region we have rather w(r) ~ 0. 
Of course the value of r,„ cannot be precisely determined for k < oo but the numerical results obtained for 
large k are quite compatible with = r^.c ~ 2.0623 and with r,„ = q. This is illustrated in Fig. [3] for 
K = 10^ where the functions w, h and the derivatives w\ h' are plotted as functions of r. The limit of the 
solutions of the equations ([5]) is therefore not differentiable at r — r^, i.e., the equations ([5]) are singularly 
perturbed about 1/k = 0. As a consequence of the absence of differentiability of the limiting solution at an 
intermediate point of the domain of integration, the numerical analysis becomes involved for increasing k. 

Coming finally to the energy of the solution, the numerical evaluation of E is fully compatible with the 
fact that the energy stays finite and that E^^ao = 2.0623. The convergence of i? for /c ^ cx) is, however, 
much slower than the convergence of the parameter q. To argue that this statement is correct, we further 
evaluate the quantity 



with the numerical profile of w(r). In fact, Ea is an estimation of the energy obtained by assuming h — 1—q/r, 
w = for r G cxd] and h ^ and the numerical value of w{r) for r G [0, g]. One should expect 
Ea{k — > oo) = Eaair^ c) since q ~ ^ in this limit. 

The value Ea is reported on Fig. [5] (see the dashed line); confirming our expectation, we find that the 
difference E — Ea tends quickly to zero (in fact exponentially fast) for k — > oo. Interestingly, Ea provides a 
reasonably good approximation to the energy E even for small values of k; for instance we find Ea/ E ~ 1.035 
for K = 1. 

5 Summary 

We have carried out a combination of analytic and numerical analysis for what we refer to as the Skyrmed 
monopole, which is a finite energy solution to the equations of the model ([1]). The asymptotic analysis 
involves the study of a one parameter family of monopole solutions, parametrised by the strength of the 
coupling of the quartic Skyrme like term. For any value of the coupling constant k we have given asymptotic 
expansions of the solutions for r ^ and r —>■ oo. We have concentrated, however, mainly on the behaviour 
of the solutions for k —> oo. The numerical analysis shows that the energy is bounded as k — s- oo, a result for 
which we have no mathematical underpinning. Given that the energy is bounded, we can however deduce 
analytically some interesting results. We find in particular that the equations ([5]) are singularly perturbed 
about 1/k = 0. This is reminiscent of Burgers' equation with a small coefficient in front of the second order 
derivative. When this coefficient goes to zero, the solutions of Burgers' equations tend to a weak solution 
with a shock, i.e., the limiting weak solution is discontinuous. In our case, the derivative of the limiting 
weak solution is discontinuous. This analytic result, and the results we get for the values of sonic typical 
constants are all supported by our numerical analysis. 
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